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SECTION  I 
INTRODUCTION 

This  is  an  interim  report  on  the  Assault  Boat  Equations,  Task  1747-04. 
Work  is  now  going  on  within  the  Naval  Training  Equipment  Center 
(NAVTRAEQUIPCEN)  to  satisfy  the  Navy's  requirement  for  an  assault  boat 
coxswain  trainer.  An  assault  boat  breadboard  device  has  been  installed  and 
is  to  be  driven  by  the  Computer  Laboratory  Sigma  7  computer  system.  The 
present  task  entails  development  and  implementation  of  a  requisite  real-time 
computer  program  based  on  landing  and  retracting  behavior  of  an  assault  boat 
from  a  beach.  Initially,  the  assault  boat  equations  and  experimental  data 
were  provided  to  the  NAVTRAEQUIPCEN  by  references  i  and  2.  The  experimental 
data,  which  were  obtained  by  Oceanics,  Incorporated  from  the  British 
Hovercraft  Corporation,  Limited,  of  England,  were  gathered  from  tank  tests 
employing  a  scale  model  of  a  typical  assault  boat  vehicle. 

The  Computer  Laboratory's  role  in  the  project  is  twofold:  first  is  the 
modification,  where  needed,  and  verification  of  the  assault  boat  mathematical 
model;  and  second  the  development  of  a  real-time  simulation  program  to  drive 
the  assault  boat  breadboard  simulator  being  built  at  the  NAVTRAEQUIPCEN. 

This  report  describes  the  in-house  work  with  the  exception  of  detailed 
software  descriptions  performed  by  the  author  under  the  direction  of  the 
project  engineer,  H.  C.  Saltzman.  Detailed  software  descriptions  will  appear 
in  internal  documentation. 
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SECTION  II 

STATEMENT  OF  THE  PROBLEM 

The  problem  was  twofold.  First,  a  real-time  simulation  program,  based 
on  the  assault  boat  mathematical  model,  had  to  be  written  to  run  on  the 
Sigma  7  computer.  Continuous  and  real-time  control  (i.e.,  instantaneous 
input/output)  of  the  assault  boat  mathematical  model  was  required  during 
verification  tests  of  the  equations  of  motion. 

The  second  problem  was  the  verification  of  the  equations  of  motion  for 
the  assault  boat.  However,  there  were  many  inconsistencies  found  between 
the  assault  boat  mathematical  model  and  the  experimental  data.  As  a  result, 
part  of  the  problem  was  to  make  modifications  to  the  assault  boat  equations 
of  motion  and  then  to  perform  verification  tests  based  on  available 
experimental  data. 


Preceding  page  blank 
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SECTION  III 
METHOD  OP  PROCEDURE 

The  procedure  was  initially  to  write  a  nonreal-time  Fortran  program 
that  ran  on  the  Sigma  7  computer.  The  purpose  of  this  was  to  make  an 
off-line  study  of  the  equations  of  motion.  Next,  a  real-time  simulation 
program  with  on-line  control  was  written.  Approximately  one-third  of  the 
program  was  written  in  Fortran  IV  and  two-thirds  in  assembly  language.  With 
no  further  addition,  the  program  has  a  cycle  time  of  60msec.  The  input  controls, 
throttle  and  rudder  angle,  were  simulated  by  the  Sigma  7  computer  rather  than 
inputted  by  the  assault  boat  hardware.  For  off-line  study,  the  card  reader 
and  line  printer  can  be  used  as  input/output  devices.  A  time  history  of 
inputs  can  be  read  into  the  simulation  program  by  the  card  reader.  Furthermore, 
the  simulation  program  (if  called  upon  by  the  operator)  will  output  up  to  eight 
variables  of  motion  to  an  analog  strip  recorder.  The  executive  program  of  the 
simulation  program  interfaces  with  the  external  sense  switches  on  the  console 
to  allow  the  operator  to  select  the  input/output  modes  described  above. 

After  the  simulation  program  was  operating  correctly,  further  tests  on 
the  validity  of  the  equations  of  motion  were  carried  out.  Because  of  the 
discrepancies  found  between  the  experimental  variables  and  the  theoretical 
numerical  values  of  motion  variables,  modifications  to  the  coefficients  were 
made,  and  in  some  cases  new  equations  were  derived  based  on  curve  fitting 
equations  to  the  experimental  data. 
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SECTION  IV 
RESULTS 

A.  SUMMARY  OF  EQUATIONS. 

At  the  tine  of  this  report,  the  following  is  a  summary  of  the  assault 
boat  equations  used.  The  glossary  contains  a  definition  of  the  variables. 


Accelerations 

^  =  m-x.  ^hull  +  Xprops  +  ^rudder  +  ^vind  +  ^waves  +  Xbeach  +  mvr^  ^ ) 

v  =  -  ^hull  +  xprops  +  ^rudder  +  xwind  +  xvaves  “  mar) 

m-y. 

Z  =  _J _  (~2jf  A*  - _©  -2jfN *  (Z+U9)  -2pg?Bz+Zyaves)  (3) 


V  =  1 


Z  =  -1  (-2j?A«33e-2iN'zz  (Z+U9) ^pg^Bz+Zyaves) 


m+2/A» 


1  =  =-L—  (K^M  I^Vves) 


Ixx-^ 


1=  1 

I„,,+2 


yy+§i3A’33 


(2jA.33(U^U2fl)4f3N.z-  |  gBi^aves> 


r  =  - 1 _  (Nhn1 1  +  NnT>nn<,  +  N. 


hull  T  “props  x  “rudder  T  “control  +  Nwind  +  Nwaves) 


Hull  Forces 

xhull  =  X1*M  +  X2VT 

(Ahead) 

(7) 

xhull  =  Xp-”*  +  x3y2 

(Astern) 

(8) 

xhull  =  +  xluv  +  x2ur 

(Ahead,  Astern) 

(9) 

Nhull  =  Nfr  +  N^ur  +  N2r|r| 

|  +  N3uv  (Ahead,  (Astern)) 

(10) 
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Propeller  Forces 

Xpropg  =  10*5no  -  .2n0u2  -  u-3  (ahead,  deep  water)  (ll) 

Xprops  =  Xpropg  +  *08un2  “  77.37u2nQ|o|j_j  (ahead,  shallow  water)  (12) 


Xprops  =  -0073u4  +  0.20u2 

(ahead) 

(13) 

NpropS  =  0.1818u4  +  5.308u2 

(ahead) 

(U) 

Xprops  =  9n02  -  n^-5  -  7.1022u2 

(astern,  deep  water) 

(15) 

X^r^ps  =  X^r^jps  .lln^u  —  62.86u^t( 

U\  (astern,  shallow  water) 

w 

(16) 

Xprops  =  6.91u  +  0 . 58u2 

(astern) 

(17) 

Vops  =  156. 91u  +  15.78u2 

(astern) 

(18) 

n0  =  mininxum(np,ng)  +  0.8|ns-np|  (astern,  ahead) 

(19) 

Rudder  Forces 

Xrudders  =  ~2.811£2 

(ahead,  astern) 

(20) 

^rudders  "  2,71u^2 

(ahead) 

(21) 

^rudders  =  0.47u2£ 

(astern) 

(22) 

Nrudders  =  -56.69u3£ 

(ahead) 

(23) 

^rudders  “  12.28u3£ 

(astern) 

(24) 

Control  Forces 

^control  =  1*25  (l(np)  —  T(ns)) 

(ahead) 

(25) 

^control  =  1*25  (T(ns)  -  T(r»p)) 

(astern) 

(26) 

T(n)  =  9.95n02 

(ahead) 

(27) 

T(n)  =  .95n03 

(28) 
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Wind  Forces 

^wind 

=  -  ^  a  Cd(A-j  Jcosxl  j+Ag  Jsin/{|)UpycosU 

(29) 

^wind 

=  Cd(A1  JcosJul+A2|sin^)URWsin>il 

(30) 

^wind 

=  0.3L  P a  CD(A1  Jccszi^lsinitjjU^. 

(31) 

(l-  |sin/ul)  sign(sirvu) 

URW 

=  /u2  +  v2 

V  RW  RW 

(32) 

URW 

=  U^cos(^y-yO-u 

(33) 

VRW 

=  Uwsin(^y-^-v 

(34) 

=  tan-1 

(35) 

Kinematics  (position  of  boat  relative  to  the  inertia  reference  frame 
at  time  tn) 


X0(tn)  =  Xc(tn_i)  +uAtcosy/  -  v^tsin  y' 

(36) 

X0(tn)  =  XQ(tn_i)  +vAtcos^+  uAtsiny' 

(37) 

At  =  tn-tn-1 

(38) 

B.  MODIFICATIONS  TO  THE  EQUATIONS  OF  MOTION  AND  TEST  RESULTS. 

The  experimental  data  used  as  a  test  criteria  were  obtained  by  Oceanics, 
Incorporated  from  the  British  Hovercraft  Corporation,  Limited.  This 
experimental  data  was  obtained  from  a  scale  model  of  landing  craft  model  (LCM) 
during  tank  tests  and  the  results  can  be  found  in  reference  1. 

Free  oscillation  tests  were  performed  to  obtain  the  natural  frequencies 
for  the  simulation  and  to  compare  these  results  with  the  experimental  values. 
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MOTION 

SIMULATION 

EXPERIMENTAL 

Heave 

2.8sec 

2.9sec 

Roll 

3.3see 

3.1sec 

Pitch 

2.1sec 

3.4sec 

Figure  1  shows  a  graph  of  these  displacements  as  a  function  of  time.  Figures 
2  through  5  show  the  variation  of  RPM  with  speed  at  self  propulsion  for 
various  conditions.  These  tests  were  carried  out  by  setting  the  RPM  of  the 
assault  boat  propellers  to  a  fixed  value  and  then  recording  the  constant 
velocity  obtained  while  maintaining  a  straight  course  with  the  rudders. 

When  these  tests  first  were  run,  the  equations  as  provided  by  Oceanics, 
Incorporated  were  erroneous.  However,  this  situation  was  rectified  by 
deriving  two  new  equations  for  the  propeller  force  in  the  X-direction  of  the 
body  axis  to  replace  the  two  equations  obtained  from  Oceanics,  Incorporated. 

xprops  =  4.93n^  +  ,82u%i0  -  .06u^  (Ahead,  deep  water) 

was  replaced  by 

Xprops  —  10 • 5n2  —  ,2u2n0  —  u? 


Also, 


Xprops  =  16.36n§  -  17.97unc  -  3.38u2 

(Astern,  deep  water) 

was  replaced  by 

Xprops  =  8. On2  —  • 9u^ • ^nQ  -  7.1022u2 

(Astern,  deep  water) 

Two  equations  for  the  propeller  force  are  required  because  the  basic 
framework  of  the  mathematical  model  treats  the  approach  to  and  retraction 
from  the  beach  as  two  separate  simulation  problems.  One  set  of  equations  is 
used  when  the  assault  boat  approaches  the  beach  and  the  other  when  the 
assault  boat  retracts  from  the  beach.  This  presents  no  problem.  However, 
Oceanics,  Incorporated  treated  propeller  force  in  the  body  X-direction  as 
if  both  the  port  and  starboard  propellers  were  rotating  at  the  same  RPM. 

As  an  outcome,  their  expression  for  combined  propeller  rotation, 
n0  =  £  (np+ns) ,  where  np  and  ng  equal  port  and  starboard  rotation  respectively, 
was  found  to  be  inappropriate  under  the  following  condition.  That  is,  when 
considering  the  net  X-body  propeller  force  due  to  a  differential  in  the  RPM 
of  the  propellers,  for  example,  just  the  starboard  propeller  rotating,  error 
in  the  forward  velocity  of  the  boat  was  introduced.  However,  another 
expression  for  propeller  rotation  was  derived  n0  =  MINIMUM  (np,ns)  +  .8|rip-ns| . 

This  expression  gave  excellent  results  as  can  be  seen  from  figures  2  through 
5.  It  will  be  noted  that  when  the  propellers  are  rotating  in  unison,  i.e., 
when  np=ns,  the  equations  for  n0  reduces  to  £  (np+ns)  which  is  the  original 

expression  provided  by  Oceanics,  Incorporated.  Furthermore,  the  new  expression 
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for  propeller  rotation  (when  used  in  the  propeller  force  equation)  gave 
excellent  results  for  the  variation  of  RPM  with  speed  when  just  the  starboard 
screw  was  rotating.  This  can  be  seen  in  figures  3  and  5. 

Figures  7  and  8  show  resistance  as  a  function  of  velocity  for  both  ahead  and 
astern1  motion.  The  original  expressions  of  resistance  to  motion  for  both 
ahead  and  astern  motion  were  inadequate.  The  inadequacy  of  the  drag  in  the 
water  became  apparent  in  deceleration  when  both  propellers  were  turned  off. 
However,  when  the  equations  of  resistance  were  altered  to  fit  the  resistance 
data  as  a  function  of  forward  velocity,  deceleration  became  more  pronounced. 

Figure  10  shows  the  variation  of  side  force  with  speed.  However,  in  order  to 
obtain  this  "good  fit",  it  was  necessary  to  derive  a  new  equation  for  the 
propeller  force  in  the  Y-body  axis  system. 

The  propeller  force  in  the  Y-body  axis  system, 

V°ps  =  7.13u+5.74u2-.75u3+.03u2, 

was  replaced  by 

^props  =  .0073u4+.20u2. 

Figure  11  shows  the  variation  of  yawing  moment  with  speed  for  ahead  motion 
and  both  propellers  rotating  in  calm  water.  As  can  be  seen,  there  is  close 
agreement  between  theoretical  and  experimental  values.  However,  as  before, 
it  was  necessary  to  derive  a  new  equations  for  yawing  moment  produced  by  the 
propellers.  The  yawing  moment  about  the  Z-body  axis  system  produced  by  the 
propellers. 


“props  =  103.13u-114.2u2+15.61u3-.75u4, 
was  replaced  by 

“props  =  .1818u4+5 , 308u2  . 

Figure  12  shows  the  variation  of  side  force  and  yawing  moment  with  speed  for 
astern  motion  and  both  propellers  rotating  in  calm  water.  The  original 
equations  provided  by  Oceanics,  Incorporated  were  not  modified  since  there 
was  close  agreement  between  theoretical  and  experimental  values.  The  dispersion 
in  experimental  data  was  due  to  aeration  effects  during  tank  tests  on  the  landing 
craft  model. 

Figure  13  shows  the  difference  between  the  port  and  starboard  RPM  as  a  function 
of  speed  when  a  straight  course  is  maintained.  However,  these  comparisons  of 
theoretical  versus  experimental  values  were  made  off-line  as  follows.  Since 
a  straight  course  is  maintained  ahead,  t  -  r  -  -fr  =  o  in  the  coupled  equations 
of  lateral  motion  (sway)  and  yaw  with  a  constant  forward  velocity,  u,  and 
lateral  velocity,  v,  maintained.  With  the  above  conditions,  the  equations  of 
motion  reduce  to  u  =  constant  (i.e.,  constant  forward  speed  and  no  yawing) 

Y-juv+Ipr0pS  =  0  (i.e.,  zero  lateral  acceleration  and  no  yawing) 
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N3urfNprops+Ncontrols  =  0  (i.e.,  zero  yaw  rate) 


where 


^controls  “  l«25(T(np)— T(ng) ), 

T(n)  =  9.95n2,  and  np,ns  =  port  and  starboard  propeller  RPS. 

By  inputting  an  initial  forward  velocity,  u,  and  starboard  propeller,  ns, 
an  equilibrium  lateral  velocity,  v,  and  starboard  port  propeller,  np,  can  be 
calculated.  This  off-line  method  was  used  to  calculate  the  theoretical  values 
for  the  difference  of  port  and  starboard  propeller  as  a  function  of  forward 
speed  when  a  straight  course  was  maintained. 

Figure  14  shows  the  variation  of  side  force  with  speed  and  rudder  angle.  The 
solid  lines  present  the  experimental  data  from  the  simulation  mathmatical 
model.  No  modifications  were  made  to  the  equation  that  provides  the  side 
force  produced  by  the  rudder. 

Figure  15  shows  the  variation  of  yawing  moment  with  speed  and  rudder  angle. 
Again,  no  modifications  to  the  equation  that  calculates  the  rudder  yawing 
moment  was  made. 
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SECTION  V 
CONCLUSIONS 

In  conclusion,  a  preliminary  simulation  assault  boat  model  has  been 
completed.  There  is  an  abundance  of  experimental  data  available  in 
reference  2  on  the  performance  of  an  assault  boat  model  in  calm  water, 
regular  waves,  and  the  landing  and  retraction  of  the  boat  from  the  beach. 
Hence,  there  is  a  good  foundation  for  improving  and  testing  the  validity  of 
further  changes  to  the  assault  boat  vehicle  equations  of  motion  in  order  to 
improve  the  simulation.  One  major  shortcoming  exists  in  the  experimental 
data;  there  are  no  acceleration  and  deceleration  curves  on  the  assault,  boat 
model.  A  major  shortcoming  of  the  math  model  is  that  although  differences 
in  port  and  starboard  propellor  rotation  are  permitted,  differences  in  the 
direction  of  propeller  rotation  is  not  permitted.  This  is  a  serious  flaw 
in  that  the  propellers  are  the  main  control  the  coxswain  has  over  the  control 
of  the  assault  boat.  This  deficiency  will  be  corrected. 
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GLOSSARY 

DESCRIPTION 

The  effective  frontal  above-water  area  of  the  boats 
hull  side. 

The  lateral  above  water  area  of  the  boat. 

Added  mass  and  damping  coefficient. 

Beam  width  of  boat. 

Average  drag  coefficient  of  the  boat. 

Acceleration  of  gravity. 

Vertical  distance  between  center  of  gravity  and 
metacenter. 

Water  depth  along  sloping  beach. 

Moment  of  inertial  about  the  x,  y,  and  z  axes. 

Moments  (i.e.,  torques)  about  the  x,  y,  and  z  axes. 
Also,  referred  to  as  the  roll,  pitch,  and  yaw  moment. 

A  roll  moment  coefficient. 

Boat  length  along  load  waterline. 

Half-length  of  boat. 

Mass  of  boat. 

Hull  force  coefficients. 

Propeller  rotational  speed,  rev/sec. 

Used  to  represent  the  effective  rotational  rate, 
rev/sec,  of  the  port  propeller. 

Rotational  rate,  rev/sec,  of  the  port  propeller. 

Rotational  rate,  rev/sec,  of  the  starboard  propeller. 

Angular  velocity  components  about  the  x,  y,  and  z  axes. 

Propeller  thrust. 

Resultant  wind  velocity  relative  to  the  boat. 

Resultant  wind  velocity  relative  to  the  inertial 
reference  frame. 
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Velocity  components  in  the  direction  of  the  x,  y  and 
z  axes. 

Weight  of  boat. 

Components  of  force  in  the  direction  of  the  x,  y,  and 
z  axes. 

Right-handed  orthogonal  system  of  moving  axes,  fixed 
in  the  body;  the  z-axis  is  directed  towards  the  bottom 
of  the  boat  with  the  x-z  axis  the  vertical  plane  of 
symmetry  of  the  boat.  < 

Right-handed  orthogonal  system  of  fixed  axes,  referred 
to  an  inertial  reference  frame  fixed  with  respect  to 
the  surface  of  the  earth.  Used  to  describe  the  position 
of  the  boat  relative  to  the  shoreline. 

Axial  hull  force  coefficients# 

Lateral  hull  force  coefficients. 

Beach  slope. 

Rudder  angle  in  radians  about  the  z-axis.  Positive 
if  rotation  is  clockwise. 

Pitch  angle,  about  the  y-axis. 

Resultant  wind  vector  angle  with  the  negative  x-axis. 

Mass  density  of  air. 

Roll  angle,  about  the  x-axis. 

Yaw  angle,  about  the  x  -axis. 

Angle  that  wind  vector  makes  with  the  inertial 
reference  frame. 


